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We study the participant-spectator matter at the energy of vanishing flow for neutron- 
rich systems. Our study reveals similar behaviour of participant-spectator for neutron-rich 
O , systems as for stable systems and also points towards nearly mass independence behaviour 

of participant-spectator matter for neutron-rich systems at the energy of vanishing flow. 
We also study the thermalization reached in the reactions of neutron-rich systems. 
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The isospin physics has attracted lot of attention of the present nuclear physics re- 
searchers around the world for the past decade. The establishment and upcoming ra- 
dioactive ion beam (RIB) facilities provide a major boon to this [1,2]. RIBs provide the 
possibility to study the nuclear matter under the extreme conditions of isospin asymme- 
try. Heavy-ion reactions induced by neutron-rich matter provide a unique opportunity 
to explore the isospin dependence of in-medium nuclear interactions, since isospin degree 
of freedom plays an important role in heavy-ion collisions through both nuclear equation 
of state (EOS) and in-medium nucleon-nucleon cross section. Role of isospin degree of 
freedom has been investigated in collective flow and its disappearance (energy at which 
flow disappears is called energy of vanishing flow (EVF) or balance energy (E^a/)) for the 
past decade. Isospin effects in flow were flrst predicted by Pak et al. [3] These isospin 
effects are due to the competition among various reaction mechanisms like Coulomb force, 
symmetry energy, isospin dependence of cross section, and surface effects. Though many 
studies are available on the energy of vanishing flow, very few studies exist in the lit- 
erature that are carried out to study other heavy-ion phenomena at the balance energy 
[4-7] . Balance energy results due to the counterbalancing of the attractive mean fleld and 
repulsive nucleon-nucleon collisions. In terms of theoretical description, it is the relative 
dominance of real and imaginary parts of G-matrix which decides the fate of a reaction. 
The dominance of nucleon-nucleon collisions at high incident energies makes the imagi- 
nary part very signiflcant. However, both real and imaginary parts of complex G-matrix 
are equally important at intermediate energies. This picture can also be looked in terms 
of participant-spectator matter and flreball concept. In Ref. [7], Sood and Puri studied 
the participant-spectator matter and nuclear dynamics for stable systems. The study 
revealed that for stable systems (N/Z ^ 1), participant-spectator matter at E^a^ is quite 
insensitive to the mass of the colliding system. It, therefore, can act as a barometer for 
the study of balance energy. No study exists in literature to demonstrate how these ob- 
servables behave for neutron-rich systems. Here we plan to extend the above study for the 
neutron-rich systems and to look whether the above participant-spectator demonstration 
still holds for systems lie far away from the stability line. The present study is carried out 
within the framework of isospin-dependent quantum molecular dynamics (IQMD) model 
[8]. 

We simulate the reactions of Ca-|-Ca, Ni-|-Ni, Zr-|-Zr, Sn-|-Sn, and Xe-|-Xe series hav- 
ing N/Z = 1.0, 1.6 and 2.0. In particular, we simulate the reactions of *^°Ca+'^°Ca 



(105), 52Ca+52Ca (85), eoCa+^oCa (73); ^SNi+^SNi (98), 72Ni+72Ni (82), ^^Ni+s^Ni (72); 
8iZr+8iZr (86), i042r+i042r (74), 1202^^1202^ (67); looSn+iooSn (82), i^ssn+i^ssn (72), 
i50Sn+i50Sn (64) and ^^^Xe+^^^Xe (76), ^^^Xe+^^^Xe (68) and ^^^Xe+^^^Xe (61) at an 
impact parameter of b/bmax = 0.2-0.4 at the incident energies equal to balance energy. 
The values in the brackets represent the balance energies for the systems. We use a soft 
equation of state along with the standard isospin- and energy- dependent cross section 
reduced by 20%, i.e. a = 0.8 o'l^'^- The reactions are followed till the transverse in-plane 
flow saturates. It is worth mentioning here that the saturation time varies with the mass 
of the system. Saturation time is about 100 (150 fm/c) in lighter (heavy) colliding nuclei 
in the present energy domain. We use the quantity " directed transverse momentum (pf^)" 
to deflne the nuclear transverse in-plane flow, which is deflned as [7, 8] 
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where y{i) andp^(i) are, respectively, the rapidity (calculated in the center of mass system) 
and the momentum of the i*^ particle. The rapidity is deflned as 

l,^|(i+M). (2) 

where E{i) and Pzii) are, respectively, the energy and longitudinal momentum of the 
i*^ particle. In this deflnition, all the rapidity bins are taken into account. It is worth 
mentioning that the E^ai has the same value for all fragments types [3]. Further the 
apparatus corrections and acceptance do not play any role in calculation of the Ehai- 

Since the balance energy represents the counterbalancing of attractive mean flled po- 
tential and repulsive nucleon-nucleon scattering and so this counterbalancing is reflected 
in participant and spectator matter as predicted in Ref. [7] In the present study, we deflne 
the participant-spectator matter in terms of nucleonic concept. All nucleons having ex- 
perienced at least one collision are counted as participant matter. The remaining matter 
is labeled as spectator matter. These deflnitions give us the possibility of analyzing the 
reaction in terms of the participant-spectator flreball model. 

In flg. 1 we display the time evolution of normalized spectator matter (left panel) and 
participant matter (right panel). The upper, middle and lower panels represent the results 
for N/Z = 1.0, 1.6, and 2.0, respectively. Lines correspond to different systems. Solid, 
dashed, dotted, dash-dotted, and short-dotted lines represent the reactions of Ca+Ca, 
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Figure 1: (Color online) The time evolution of spectator matter (left panels) and partici- 
pant matter (right panels) for systems having N/Z = 1.0, 1.6 and 2.0. Lines are explained 
in the text. 



Ni+Ni, Zr+Zr, Sn+Sn, and Xe+Xe, respectively. From figure, we find that at the start 
of the reaction, all nucleons constitute the spectator matter. Therefore, no participant 
matter exists at t = fm/c. As the reaction proceeds we have the decrease in spectator 
matter with corresponding increase in participant matter. We also find that for lighter 
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Figure 2: (Color online) The N/Z dependence of participant and spectator matter. Sym- 
bols are explained in the text. 



systems like Ca+Ca and Ni+Ni, the transition from spectator to participant matter is 
swift and sudden whereas for the heavier colliding nuclei, the transition is slow and gradual 
as predicted in Ref. [7]. This is because of the fact that lighter reactions occur at relatively 
high energies. At the end of the reactions, we have nearly the same participant matter, 
which indicates the universality in balancing the attractive and repulsive forces. We 
also see that similar behaviour exists for all N/Z ratios. This indicates that participant- 
spectator behaviour is similar for neutron-rich systems as for systems lying on the stability 
line (N/Z = 1). 

In fig. 2 we display the N/Z dependence of participant and spectator matter. Upper 
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Figure 3: (Color online) The system size dependence of participant and spectator matter 
for different N/Z ratios. Various symbols are explained in the text. 



(lower) panel displays the spectator (participant) matter. Squares, circles, triangles, di- 
amonds, and pentagons represent the reactions of Ca+Ca, Ni+Ni, Zr+Zr, Sn+Sn, and 
Xe+Xe, respectively. We find that for all the system masses participant-spectator mat- 
ter is almost independent of N/Z. There is a very slight increase (decrease) in spectator 
(participant) matter with N/Z of the system. 

In fig. 3, we display the system size dependence of the participant and spectator 
matter. Open (solid) symbols represent participant (spectator) matter. Upper, middle 
and lower panels represent the results for N/Z = 1.0, 1.6 and 2.0, respectively. We see 
that participant-spectator matter follows a power law behaviour (oc A^) with the system 
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Figure 4: (Color online) The time evolution of the anisotropy ratio (upper panel) and 
relative momentum (lower panel) for various systems having N/Z = 1.0. Lines have the 
same meaning as in Fig. 1 



mass. The power law factor is -0.12± 0.04 (0.02± 0.01), -0.14± 0.06 (0.02± 0.01), and 
-0.03± 0.06 (0.01± 0.01) for spectator (participant) matter having N/Z = 1.0, 1.6 and 
2.0, respectively. Thus, a nearly mass independent behaviour is obeyed by the participant 
and spectator matter for all the N/Z ratios. 

In fig. 4 we display the time evolution of anisotropy ratio < Ra > (upper panel) and 
relative momentum < Kji > (lower panel) for different system masses having N/Z = 1.0. 
The < -Ra > is defined as 

VPx + \/Pv , s 

< Ra >= ^ " ^ ' . (3) 



This anisotropy ratio is an indicator of the global equilibrium of the system. This 
represents the equilibrium of the whole system and does not depend on the local positions. 
The full global equilibrium averaged over large number of events will correspond to < /?„ > 
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Figure 5: (Color online) The system size dependence of anisotropy ratio for various N/Z 
ratios. 



= 1. The second quantity, the relative momentum < K^ > of two colliding Fermi spheres, 
is defined as 



<Kr> = < \Pp{f,t)-PT{f,t)\/h>, 



(4) 



where 



P^ir,t) 



j:f=,mpAr,t) 



i = l,2. 



(5) 



Pj{f,t) 

— * 

Here Pj and pj are the momentum and density of the j'th particle and i stands for either 
projectile or target. The <Kji> is an indicator of the local equilibrium because it depends 
also on the local position r. 



From figure 4(a) (upper panel), we see that anisotropy ratio increases as the reaction 
proceeds and finally saturates after the high density phase is over. We also see that the 
influence of system size is very less on anisotropy ratio and hence indicates towards the 
equilibrium of the system. From fig. 4(b) (lower panel) we see that relative momentum 
decreases as the reaction proceeds and the smaller value of <Kf>> at the end of the 
reaction indicates toward the better thermalization of the matter. We also see from the 
figure that <Ra> ratio saturates as soon as high density phase is over which signifies 
that the nucleon-nucleon collisions happening after high density phase do not change the 
momentum space significantly. 

In fig. 5 we display the system size dependence of anisotropy ratio for systems having 
N/Z ratios 1.0, 1.6 and 2.0. From figure we see that anisotropy ratio follows a power law 
behaviour (oc A"^) with system size. The power law factor is 0.07 ±0.03, 0.04 ±0.008, and 
0.07 ±0.009 for N/Z ratios 1.0, 1.6 and 2.0, respectively. 

In summary, we studied the participant-spectator matter at the energy of vanishing 
flow for neutron-rich systems. The study revealed a similar behaviour of participant- 
spectator for neutron-rich systems as for stable systems and also pointed towards a nearly 
mass dependence behaviour of participant-spectator matter of neutron-rich systems at 
the energy of vanishing flow. Similar mass independent behaviour is also found for the 
anisotropy ratio. 

This work has been supported by a grant from Centre of Scientiflc and Industrial 
Research (CSIR), Govt, of India and Indo- French Centre For The Promotion Of Advanced 
Research (IFCPAR) under project no. 4104-1. 

References 

[1] W. Zhan et al, Int. J. Mod. Phys. E 15, 1941 (2006); see, e.g. 
http://www.impcas.ac.cn/zhuye/en/htm/247.htm. 

[2] Y. Yano, Nucl. Inst. Methods B 261, 1009 (2007). 

[3] R. Pak et al, Phys. Rev. Lett. 78, 1022 (1997); ibid. 78, 1026 (1997). 

[4] Z. Y. He et al., Nucl. Phys. A598, 248 (1996). 

[5] B. A. Li et al, Phys. Rev. C 48, 2415 (1993). 



[6] S. Soff et al, Phys. Rev. C 51, 3320 (1995). 

[7] A. D. Sood and R. K. Puri, Phys. Rev. C 70, 034611 (2004); tbid. C79, 06461^ 
(2009); tbid. C73, 067602 (2006); tbid. Int. J. Mod. Phys. E 15, 899 (2006). 

[8] C. Hartnack et al, Eur. Phys. J. A 1, 151 (1998). 



10 



